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Introduction 

£ 

There are different means which permit determination of /24 

the form and extent of ice accretion on aircraft airframes and of 
qualifying deicing and antiicing systems : 
natural icing flight tests; 
artificial icing flight tests; 

of actual parts installed on a carrier aircraft behind 
a water atomizing grid; 

of the aircraft itself, flying in an artificial cloud 
behind a tanker aircraft; 
fixed point ground tests (helicopters); 
wind tunnel tests; 

actual parts or full scale parts; 
complete models or reduced scale models. 

These latter wind tunnel test techniques have been used by 
ONERA (National Office of Aerospace Studies and Research) since 
1962, in the large, naturally cooled Modane SI wind tunnel [1], 
during the severe winter months. The first test, carried out on 
a full scale part of the leading edge of a Vickers-Arms trong VC 10 


Numbers in the margin indicate pagination in the foreign text. 


1 



OFp^ PAGB © 
P0 °* QUALIIV 


wing, was followed in later years by other tests on full scale 
parts of various aircraft : 

wing and stabilizer parts of a Nord 262 (now known as the 
Aerospatiale Fregate) , of a Potez 840 and of a Breguet 941 S STOL; 
radome of a C 160 cargo Transall. 


To perform these tests, test section No. 1 of the large SI 
wind tunnel was equipped with a grid, which permits production of 
a homogeneous ice fog ona0.8mx2m area. 

This atomizing grid was supplied by Vickers-Armstrong Ltd., 
which previously used it on its Vanguard flight test bed. The 
90 Napier and Sons, Ltd. injectors, supplied with warm water 
and compressed air, provide a cloud of supercooled water droplets, 
the median diameter of which can be regulated between 15 and 25 yin, 
and a liquid water content between 0.2 and 2.5 g^m - ^ f or speeds 
up to 100 m*sec~l' an( ^ with some limitations, up to 120 m*sec“l. 

The problem of icing in the wind tunnel has been considerably 
modified by the study of ice accretion on the Concorde supersonic 
transport in the holding, approach and landing flight envelopes 
under icing conditions. In fact, this aircraft has an ogive delta 
wing, the vortex type flow of the top surface of which is closely 
connected to the shape of the entire wing and the presence of the 
front portion of the fuselage. Therefore, in order to determine 
the areas, shapes and extent of icing, it is necessary to perform 
tests on a complete airframe. Given the dimensions of the air- 
craft, full scale tests are possible only in actual flight. Wind 
tunnel icing tests performed at NRC Ottawa on models or semimodels 
of flat delta wings [2, 3], have demonstrated large ice accretions 
on the top of the wing, and they permit the belief that there is 
a prohibitive accretion of weight for the aircraft, as well as 
perturbation of the longitudinal stability. A collecting study 
in the hydrodynamic tunnel, carried out at the University of 
Bristol from 1963 to 1967 [4], led to the same qualitative conclu- 
sions. However, these models were not representative of the Con- 
corde aircraft, with its cambered and twisted ogive fuselage and 
wing section. Therefore, it was necessary to decide on a larger 
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scale. In order to have a better determination of the shapes and 
extent of the ice and the risk of perturbations caused by the ice 
accretion, so as to intervene as quickly as possible, at the level 
of aircraft fabrication, in the situation, extent and power of 
deicing systems. The Official English and French Services, in 
liaison with the British Aircraft Corporation (BAC) , which made a 
detailed study of the phenomenon, decided to perform 1/6 similitude 
icing tests on a half model of the supersonic transport, in wind 
tunnel SI of the ONERA Center at Modane-Avrieux . ONERA was respons- 
ible for definition and production of the appropiate test 
installation. 

The 1/6 similitude icing was extended later to 1/12 similitude, 
for application to other tests. ORIGINAL PAGE 

OP POOR QUALITY 

Theoretical Study of Laws Governing Similitude Icing 

Theoretical study of the laws governing similitude icing 
requires analysis of various physical phenomena. The first of 
these phenomena is the mechanical collection of water drops, 
which brings into play research on the trajectories of these 
drops in the aerodynamic field of the flow surrounding the ob- 
stacle being iced. A study of this problem was accomplished as 
early as 1946, by Irving Langmuir and Katherine B. Blodgett of 
the General Eleectric Company [53. 

The second of these phenomena is the distribution of water 
and, then, ice on the profiles. This brings into play the concen- 
tration 'of water of the clouds, the icing time and, also, the 
convective heat transfer coefficient, which was defined by Thomas 
F. Gelder and James P. Lewis in [6]. The third phenomenon and, 
probably, the most difficult to deal with, is thermal equilibrium 
of a wall under icing conditions. A very detailed study of this 
problem was performed by Bernard L. Mes singer of the Lockheed 
Aircraft Corporation [7], in which he brought the concepts of /2$ 

the relative heat factor b and the freezing fraction n. 
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These various studies were joined together, so as to learn 
the relations between the parameters which govern similitude icing 
in the wind tunnel, by Hauger, K. G. Englar and W. W. Reas'er of 
Douglas Aircraft Company, Inc. [8]. 

For the similitude icing tests on the Concorde, as well as 
for the similitude icing tests on simple models intended to con- 
firm theoretical calculations, the parameters of the icing clouds 
were calculated, beginning with the work done by R. Googan, E. T. 
Jackson and J. H. Hubbold of the British Aircraft Corporation [9, 
10]. This work permitted selection of the wind tunnel in which 
the tests should be performed, and accomplishment of the trans- 
position of the holding and approach flight envelopes of the 
aircraft to that being simulated in the wind tunnel, within the 
framework of limited icing similitude. This similitude is limited 9 
because it is not possible to observe simultaneously the Mach 
number and Reynolds number conditions of the aerodynamic flow 
around the obstacle and the Reynolds number conditions for the 
droplets. However, it is found that, under flight conditions in 
which icing can be a hazard, the aerodynamic field in the collecting 
zones is not sensitive to these two parameters. 

Analysis of the aerodynamic and inertial forces on the water 
droplets- in the vicinity of bhe obstacle permits formula (1) to 
be obtained, which connects the various parameters of the similitude 
of the droplet trajectory 


a m \ 2 X / Voo m \ 

a G / ' Voo(3 J 
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1 (-r.G + m) 


In this formula, x is the exponent of the Reynolds number in 
the function which defines the ratio between the actual drag of a 
droplet and that given by Stokes law. Fig. 1 shows that, for a 
droplet Reynolds number area between 6 and 120, which also covers 
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the area of flight of the aircraft, as well as 1/6 and 1/12 scale 
tests, the value of x represented by the slope of the best straight 
line through the points can be assumed to be 0.39. In most cases, 
ratio T am /T a Q is about 1, and the practical forumla to use for the 
calculations becomes 
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( 2 ) 


Analysis of the various parameters acting on the volume and 
distribution of icing on the profiles, as well as introduction of 
the concepts of the fraction of water changing to the ice state 
n and of the relative heat factor b [7] permits formulas (3) and 
( 4 ) to be obtained. They define, respectively, the ratio of the 
liquid water content of the fog (LWC) and the icing time ratio (x) 



The hypotheses on which these equations are based, i.e., that 
the relative heat factors and the icing fractions are equal in 
flight and on the model (b m =bQ and n m =nQ) , can be satisfied or 
not, according to the values of the flow variables. Therefore, 
it is advisible to determine the relationships between the vari- 
ables which ensure equilibrium between the convective, evaporative 
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Fig. 1. Ratio between actual drag co- 
efficient of droplet and that given by 
Stokes law. 


and kinetic heat quantities. 

The formulas obtained in [7] 
give the equation which con- 
nects the air temperature, 
speed and pressure to the two 
parameters n and b, in the case 
of an unheated wall 

1 0,58 x 10 !i _ (5) 

Pa 


P Voo 

la (I f b) + 1732 + 79,7 nb + (3,6458 + b)— ~-r- 

Pa 8J7 j 


If a point of the aircraft flight envelope is defined by 
p aG 5 t aG> P ooG and Voo G» relationship between n and b is defined. 

Then, mathematically, it is sufficient, with the wind tunnel /26 

relationships known, which connect P am and V«, m , as well as t am and 
Poo m , to satisfy equation (5) so that the relationship between n and 
b is the same as that in the case of flight. Satisfaction of this 
relationship determines a pair of values of V°° m and t am which, 
substituted in equations (2), (3) and (4), permits determination 
of the homologous values of all the similitude icing parameters of 
the point of flight under consideration. 


Test Resources 

The specifications necessary to the British Aircraft Corpora- 
tion to satisfy simulation of the holding and approach flight 
envelopes of the Concorde aircraft under icing conditions [10] 
required completely redoing the ice fog installation in test 
section No. 1 of wind tunnel SI. These are the limiting conditions 
of these specifications, which defined the performance of the in- 
stallation. 

The new atomizing grid can hold 444 air fractionation injec- 
tors, distributing a homogeneous cloud of droplets, the median 
diameter of which can be controlled between 10 and 20 ym on a 
2 m x 2 m area (see Fig. 2). The liquid water content is 
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adjustable between 0.4 and 10 
g*m“3, for speeds which can 
reach 100 nr sec" 1 . 

The air and water sup- 
plies of this grid can be 
controlled automatically, 
to perform alternating icing 
(maximum discontinuity condi- 
tions and maximum continuous 
or clear sky conditions alter- 
nately, according to official 
specifications), while holding 
the droplet diameter practi- 
cally constant, with periods 
of f ’om 3 seconds to several 
tens of seconds. Study of 
the aerodynamic flow behind 
the grid was carried out at 
reduced scale by B.A.C., in 
Filton wind tunnel No. 4 and, 
then at full scale in wind 
tunnel SI by ONERA . The loss 
of positive or negative charge 
according to the atomizing air 
flow, was determined down- 
stream of the grid. 


Fig. 2. Icing grid and G14 model of 
SST installed in wind tunnel SI MA. 


The preliminary tests 
demonstrated that the water 
droplets of the fog evaporate if the air in the flow is not sat- 
urated. This evaporation changes the droplet diameter and the 
liquid water content of the ice fog, which then produces ice 
accretions which are not representative of the assigned conditions. 
To alleviate this disadvantage, wind tunnel SI was equipped with a 
saturation unit, consisting of eight atomizing tubes, which can 
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Fig. 3. Icing and saturation installation 
of return wind tunnel SI, Modane-Avrieux 
Center. 

Key: a cold water 

b grid 
c hygrometer 
d warm water 
e saturated 
f water temperature 
g air temperature 
h telemetry cabinet 
i air pressure 
j drying air 

k three way electric valves 
1 electric valves 
m coarse and fine flow regulators 
n fine flowmeter 
o compressed air 
p air reheater 
q coarse flowmeter 
r mixer 


discharge 900 Jl*hour _1 each, 
thus permitting saturation 
in the test section by 
raising the relative humid- 
ity from 60 to nearly 100$. 
The relative humidity is 
monitored continuously, 
with an automatic contin- 
uous sampling dew point hygro 
meter in the test section. 

The 8 m diameter test 
section of wind tunnel SI 
permits the testing of large 
models with a small ob- 
struction. Access to these 
models is very easy five 
minutes after icing is 
stopped. This easy access 
permits recording of the 
extent and shape of the ice, 
by photography at different 
distances and by direct 
measurement on the ice 
accretions. Interpretation 
of the photographs is facili- 
tated by the use of cross 
ruled sights or by very 


distinct reference marks painted on the models, to mark the chord 


percentages or angles. The test section floor has a balance, 
which permits setting the angle of attack and measurement of the 
aerodynamic forces of the models before and after icing, as well 


as a compressed air inlet to supply the jet with primary air, 
which permits simulation of the air inlet flow coefficients of 


some models. 
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The test section walls can be equipped with projectors and 
high speed cameras (1200 images • sec -1 ) , required for study of the 
trajectories of ice fragments detached naturally or artificially 
from the upstream parts of the models. 



Similitude Icing Tests ' 

Still with no element of comparison between inflight icing 
and icing in the wind tunnel at reduced scale, it was necessary 
to support the similitude icing calculations by full scale and 
reduced scale comparative icing tests in the wind tunnel. There- 
fore, ONERA, in parallel with the manufacturer icing tests, carried 
out series of similitude icing tests, limiting the study, either 
to simple shapes such as cylinders, or to wing section shapes 
close to those used in the test for which these experiments were 
conducted. The tests first were conducted within the framework 
of 1/6 similitude and, then, extended to 1/12 similitude. In 
this case, the calculations were the same as for the similitude 
between flight and the wind tunnel, but conduct of the two tests 
in one wind tunnel, in which neither the pressure nor the tempera- 
ture can be regulated, does not permit rigorous achievement of 
the theoretical thermodynamic similitude of icing. 

Each similitude icing point involved two tests, one full 
scale (subscript G) and the other at scale k (subscript m) . The 
mean droplet diameters used was &q=20 ym (standard /value given by 
the specifications) and a m =10 ym (minimum certain value permitted 
by the installation) , or a a m /aQ ratio =1/2, similar to those used 
in the flight-wind tunnel similitude. 

Here, for example, are the characteristics of two similitude 
icing points, one in 1/6 scale and the other, 1/12, carried out 
during the tests. 

1/6 Similitude Icing Test 

Full Scale Icing : 

Vq= 45 m'sec™ 1 , t a Q=-5.4° C, aQ-20 ym#, L.W.C.q= 4 g*m"^, tq=10 min. 
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1/6 Scale Icing: 


V m =15 ni'sec" 1 , t am = -5 0 C, a m =10 yml, L.W.C. m =7.2 T m =2.75 min. 

1/12 Similitude Icing Test 

Full Scale Icing : 

Vq= 76.2 m-see -1 , C, aQ=20 ym#, L.W.C.q= 2.4 g.*m'*3 J Tq = 

10 min. 

1/12 Scale Icing : 

V ra =8 m’sec -1 , t am =-5° C, a m =10 pm#, L.W.C. m =6 g*m~3 j x m =3 min. 

It is very difficult to determine angles or dimensions on 
the models with a precision better than 10 to 15 % • In the case 
of full scale icing, this is due to the branched forms of the 
deposits and, in the case of small scale icing, to the small 
actual thickness of the deposits. 

In the case of icing where the similitude rules were not 
observed, very dissimilar ice deposits appear on the models, 
both from the point of view of extent and shape of this ice. 

Fig. 4 shows a set of photographs, in which the photographs on 
the left a. represent icing obtained on models associated with 
accurate clouds and those on the right b. icing obtained with 
models and clouds which do not observe the similitude rules. 

The icing on the cylinders of 4.1. a and .b exhibit very dif- 
ferent icing angles (80 and 140°). As to the thickness of the 
ice reduced to the same scale (full scale), it differs very much 
(30 mm and l8o mm). On the delta wing sections of 4. 2. a and .b, 
the icing on the leading edge, which develops in the case of a 
normal test, scarcely appears on the large model. In the case of 
icing on wing sections of thick cross section, in 4. 3. a 
and .b, the differences are still greater: proper icing produces 
ice extending up to 2.5 % of the chord, while the ice on the small 
model extends up to 15 % of the chord, with very different accretion 
forms . 
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Fig. 5. 1/6 similitude icing on cylinders, 


ORIGINAL PAGFJ^ 
0F POOR QUALlTV 

Figures 5 and 6 show 
the results of 1/6 and 
1/12 similitude Icing, in 
the form of photographs 
and directional diagrams, 
obtained on cylindrical 
specimens . 

The icing angles of 
the ice accretions are 
practically indentical: 

65 to 75° for 70° in the 
case of 1/6 icing and 55 
to 65° for 60° in the case 
of 1/12 icing. The ice 
thickness reduced to the 
same scale is similar: 
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The small scale icing is more regular than that at full 
scale, where lateral projections appear. On the other hand, in 
all cases, the localized area around the stopping point of the 
flow is formed of transparent ice. 

Figure 7 shows the shape and extent of ice obtained at 1/6 
scale, on two 75° swept back delta half wings, set at a 15° angle 
of attack. Ice only develops on these wings on the bottom side 






Fig. 6. 1/12 similitude icing on cylinders. 



Fig. 7. 1/6 similitude icing of delta 

half wings. 
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of the leading edge, in 
the form of a laminar 
accretion, of which a 
section perpendicular to 
the leading edge has the 
shape of an isosceles 
triangle attached by one 
of its apexes. This 
accretion shape also has 
been observed at 1/6 scale, 
as well as at full scale. 

The Figure 7 photo- 
graphs show a relatively 
larger accretion on the 
full scale model, than 
on the 1/6 scale model. 

The ice thickness is 55 mm 
and 6.5 x 6=39 mm, re- 
spectively . 

Figure 8 shows the 
mounting of two models of 
full scale and 1/12 scale 
thick cross section wing 
sections in the test sec- 
tion. The right portion 
of this figure shows a 
photographic montage, 
which has reduced the two 
models to the same size, 
guided by the spacing of 
the 20, 25, 10, 5 and 2.5 
chord percentag markers. 
The similarity of the 
accretions is noteworthy, 
both as to thickness 
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(23 mm and 1.7 x 12=20.4 
mm) , as well as for the 
extent of the Ice (2.5% 
of chord). This test 
was performed at 0 angle 
of attack. The icing 
shown in Fig. 9 on the 
same models, was carried 
out at an 8° angle of 
attack. On the top views 
(left photos), no ice 
accretion appears down- 
stream of the 2.5% chord. 
On the bottom views 

(right photos), small accretions appear far downstream of the 
leading edge, and they only become large on the two models between 
5 and 2.5%. The two tests again demonstrate good consistency. 










Fig. 8. 1/12 similitude icing on thick cross 

section wing sections at 0 angle of attack. 


Altogether, these tests, performed under conditions observing 
the icing similitude rules, demonstrate the good agreement obtained, 
both as to the extent of the accretions and their shapes. Never- 
theless, in general, a slight lack of ice must be noted in the 
case of small scale icing, which perhaps can be explained by the 
production of a certain percent of ice crystals, at the level of 
atomization of the water of the fog. This percentage of crystals 
is larger in small scale tests. This requires finer atomization 
and higher air pressure. 


The tests also have sh own that, in tests performed at con- 
stant temperature, it is not necessary to reproduce the internal 
structure of the aircraft, in order to obtain the same ice accre- 
tions, and that the material of which the models are made has no 
appreciable effect on the accretions (duralumin, steel, araldite, 
wood). This is very important, for it then is possible to per- 
form reduced scale icing tests on the classical models used for 
aerodynamic tests. Most certainly, these tests do not involve 
the requirement of determination of the efficiency of equipment 
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Fl§* 9. 1/12 similitude icing on two thick cross section wing 

sections at 8° angle of attack. 

such as de-icers or antiicers, but, rather, of sticking to deter- 
mination of the areas, extent and shapes of the ice. 

The effectiveness of saturation of the test section has been 
verified, by performing icing tests with a variable distance 
between the atomizing grid and the models. There was no difference 
in the range of distances tested (from 4 to 10 m) . These tests 
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not only show that there is no evaporation of the water droplets, 
but also that the parameters of the ice fog are sufficiently 
accurate and constant, that no difference develops between tests 
performed under the same rated conditions. 

Conclusions 

The National Office of Aerospace Studies and Research has, 
in the large SI wind tunnel of its Modane-Avrieux Center, equipped 
with its new ice fog producing installation, one means of investi- 
gation necessary to builders, to specify the icing hazard which 
can be encountered in flight and to define protection systems. 

In recent years, ONERA, in close collaboration with the 
British Aircraft Corporation, the National Industrial Aerospace 
Company and the official services, has extended the field of 
use of the SI wind tunnel to icing, by successfully studying and 
testing the similitude icing technique. Therefore, it is now 
possible to determine, long before producing the aircraft and on 
the models normally used for classical aerodynamic tests, the 
areas, shapes and thickness of ice, before it is formed on an 
aircraft in flight under icing conditions. 

Some limitations, such as the annual useable cold period, /31 

temperature minimums and limited similitude, means that this 
test resource cannot be completely substituted for all the others. 
Nevertheless, the flight test program of an aircraft can be consid- 
erably eased, if the wind tunnel tests cut across typical points 
of this program. In particular, the most critical cases of icing, 
which are dangerous in actual flight, can be simulated in a wind 
tunnel without any risk, and their results can be extrapolated to 
flight, based on comparisons in less dangerous cases. 
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